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ABSTRACT 

We propose a possible new model for the formation of chain galaxies at high redshift. Our model is 
summarized as follows. Step I: Successive merging of subgalactic gas clumps results in the formation of 
a galaxy with a mass of 1O 11_12 M0 at redshift z ~ 5. Step II: Subsequently, supernova explosions occur 
inside the galaxy and then blow out as a galactic wind (or a superwind). This wind expands into the 
intergalactic space and then causes a large-scale shell with a radius of several hundreds of kpc. Since 
this radius may be smaller than the typical separation between galaxies, interactions of shells may also 
occur, resulting in the formation of a large-scale gaseous slab. Step III: Since the shell or the slab can be 
regarded as a gaseous sheet, filament-like gravitational instability is expected to occur. Step IV: Further 
gravitational instability occurs in each filament, leading to intense star formation along the filament. 
This is the chain galaxy phase. Step V: The filament collapses gravitationally into one spheroidal system 
like an elliptical galaxy within one dynamical timescale of the filament (~ 10 8 yr). Therefore, it seems 
quite difficult to find remnants of chain galaxies. 

We also discuss that shocked shells driven by superwinds may be responsible for some Lyman limit 
systems and damped Lya systems because their H I column densities are expected to be A^ji J> 10 19 
cm . 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: starburst - stars: formation - 
quasars: absorption lines - intergalactic medium 



1. INTRODUCTION 

Recent high-resolution optical imaging studies of some 
deep fields have revealed a new population of galaxies at 
high redshift; chain galaxies [Cowie et al. 1995 (hereafter 
C95); van den Bergh et al. 1996 (hereafter vdB96)]. C95 
found 28 candidates of such chain galaxies in their deep 
HST I band (F814W) images of the two Hawaii deep sur- 
vey fields, SSA 13 and SSA 22. vdB96 also identified such 
chain galaxies in the Hubble Deep Field (HDF: Williams 
et al. 1996); HDF 2-234 (the tadpole or head-tail galaxy) 
and HDF 3-531 (the chain galaxy). The chain galax- 
ies tend to be very straight in morphology; their average 
major-to-minor axial ratio is ~ 5. Another interesting 
property of the chain galaxies is that they are bluer on 
average than galaxies with similar / magnitudes. Based 
on these observational properties, C95 proposed that the 
chain galaxies comprise a new population of forming galax- 
ies at high redshift [z ~ 0.5 — 3) because there is no local 
counterpart. 

However, Dalcanton & Shectman (1996) proposed a new 
interpretation that the chain galaxies are edge-on low sur- 
face brightness (LSB) galaxies at low redshift (z ~ 0.03). 
They first estimated true axial ratios of the chain galaxies 
studied by C95. They showed that the true axial ratios 
range between 2 and 14 in most cases but some extreme 
cases have the ratios higher than 20 if the HST images 
are deconvolved with the typical point spread function. 
These higher axial ratios may favor the edge-on disk inter- 
pretation. They also showed that the ellipticity distribu- 
tions of low-redshift LSB galaxies are consistent with the 
low-rcdshift chain galaxies being edge-on manifestations 
of disks using their database of Las Campanas redshift 



survey LSB galaxies. 

Recently, Bunker et al. (2000) obtained Keck opti- 
cal spectroscopy of the Hot Dog galaxy [= HDF 4-555.1 
(Williams et al. 1996) = C 4-06 (Steidel et al. 1996)] at 
z — 2.803, which is one of chain galaxies in the HDF. They 
found that any systemic rotation of this galaxy is fairly 
small; < 100 km s _1 . This suggests strongly that this 
galaxy is not an edge-on disk galaxy. Further, the pres- 
ence of N v and He n emission lines needs massive O stars 
as their ionization sources at least. Its blue color also indi- 
cates that the stellar populations should be young (< 10 s 
yr). This study reinforces the possibility that some chain 
galaxies are really forming galaxies with a linear structure 
at high redshift as suggested by C95 (see also Abraham 
1997). Even if the majority of chain galaxies might be 
edge-on LSB galaxies at low redshift, it seems urgent to 
investigate possible origins of chain galaxies at high red- 
shift such as the Hot Dog galaxy. 

2. OBSERVATIONAL PROPERTIES OF CHAIN GALAXIES 

Prior to presenting our new model, we summarize im- 
portant observational properties of chain galaxies found 
by C95 and vdB96. A brief summary of the two surveys 
is given in Table 1. 

2.1. The C95 Sample 

C95 found 28 candidates of chain galaxies in the two 
Hawaii deep survey fields, SSA 13 and SSA 22. The 
chain galaxies tend to be very straight and their major-to- 
minor axial ratios range from 1.8 to 9.5 with an average 
of 4.8 ± 2.0 for 21 galaxies brighter than / = 25. Two 
of them are intermediate-redshift galaxies; z — 0.489 for 
Chain and z — 0.505 for Chain 1. Another two galaxies 
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are located at z = 1.36 (SSA 22-16) and at z ~ 2.4 (SSA 
22-24) although the latter redshift estimate is based on 
interpreting a break in the continuum at 4000 A as the 
onset of intergalactic Lya forest. 

It is found that the chain galaxies are bluer on average 
than galaxies with similar / magnitudes studied by them. 
In particular, they are very blue in B — K , suggesting that 
they are not relatively normal galaxies in the rest ultra- 
violet band and that the peculiar morphologies are not a 
consequence of the distribution of the star- forming regions. 

2.2. The vdB96 Sample 

vdB96 classified eleven galaxies either as the chain or 
as the tadpole galaxies. Their basic data are given in Ta- 
ble 2. In order to estimate their linear sizes and absolute 
magnitudes, we need their redshifts. However, a spectro- 
scopic redshift is available only for the galaxy 3-379, the 
ID number in vdB96 (Steidel et al. 1996): z = 2.775. 
It is interesting to note that this redshift is very close to 
that of the Hot Dog galaxy, 2.803 (Steidel et al. 1996; 
Bunker et al. 2000). For galaxies in the HDF, Fernandez- 
Soto, Lanzetta, & Yahil (1999) estimated their photomet- 
ric redshifts by using the SED (spectral energy distribu- 
tion) fitting method 1 . We use these photometric redshifts 
for the remaining 10 galaxies. Unfortunately, photometric 
redshifts are not reliably determined for the two galaxies; 
3-367 and 4-627, which are not considered further. Then 
we obtain a sample of the nine chain galaxies. We use this 
sample to investigate the observational properties of chain 
galaxies; although the chain galaxies studied by C95 con- 
stitute the largest sample, their photometric data are less 
available than those found in HDF. 

The redshifts of the nine galaxies ranges from 0.5 to 2.8, 
being similar to the range estimated for the C95 sample 
of chain galaxies. Adopting a cosmology with the Hubble 
constant, H = 100 h km s _1 Mpc -1 where h — 1 and the 
density parameter £! = 0.01, we estimate the linear size 
(d) and the absolute V magnitude (My) for each galaxy. 
The results are given in Table 2. We find that d ranges 
from 1.2 kpc to 3.8 kpc and My ranges from —16.2 to 
—21.5. Their average values are < d >~ 2.4 ± 0.8 kpc 
and < My >~ —19.0 ± 1.7, respectively. Since any incli- 
nation correction has not been made for the linear sizes, 
they are lower limits. If they are disks viewed edge-on, the 
linear sizes appear smaller than those of the present-day 
galaxies. On the other hand, the absolute V magnitudes 
are almost comparable to those of the present-day galax- 
ies. Frequency distributions of the redshifts, sizes, and 
absolute magnitudes are shown in Figure 1. 

The major-to-minor axial ratios of the nine galaxies 
range from 1.6 to 5.9 and their average is 3.1 ± 1.3, be- 
ing similar to those of the chain galaxies studied by C95. 
However, as claimed by Dalcanton & Shectman (1996), 
these axial ratios should be typically larger than 10 if the 
images are deconvolved with the point spread function. 

In summary, the vdB96 chain galaxies have the fol- 
lowing important observational properties; d ~ 2.4 kpc, 
My 19 (or L v ~ 3.4 x 10 9 i Q ), and z ~ 0.5 - 3, all 



of which are similar to the chain galaxies of C95. These 
properties reinforce that a population of chain galaxies 
does actually exist at intermediate and high redshift, even 
though edge-on LSBs at low redshifts may share this ap- 
pearance. 

3. MODEL 
3.1. Background and Motivation 

A simple explanation for high- z chain galaxies may be 
that linear structures form during the collapse of the pro- 
togalactic gas (C95; Bunker et al. 2000). C95 suggested 
that when star formation turns on it triggers induced star 
formation along the line of maximum density, being anal- 
ogous to sequential propagation of OB star associations 
in molecular clouds. On the other hand, Bunker et al. 
(2000) suggested that star formation could be triggered 
by collapse along a filament (see also Abraham 1997). An 
alternative idea may be string-induced wakes as noted by 
C95. In this paper, we propose another idea that chain 
galaxies arise from gaseous filaments which are formed in 
a large-scale shocked shell caused either by a superwind or 
by interaction between superwinds of forming galaxies. 

The starting point of our model is that the superwind 
activity occurs in some giant galaxies which were assem- 
bled from a number of gaseous clumps at high redshift 
within the framework of so-called cold dark matter scenar- 
ios. This assumption implies that initial starbursts occur 
in some giant galaxies in their forming phase. We note 
below that a couple of lines of evidence for such luminous 
forming galaxies have been obtained recently. 

(1) A number of faint dusty galaxies have been found by 
deep surveys at submillimctcr (Smail et al. 1997, 
1998, 1999, 2000; Hughes et al. 1998; Barger et 
al. 1998, 1999a, 1999b; Barger, Cowie, & Sanders 
1999; Barger, Cowie, & Richards 2000; Ivison et al. 
2000b) and at far infrared (Kawara et al. 1998; 
Eales et al. 1999; Puget et al. 1999). Most of 
these objects are considered to be located at inter- 
mediate and high redshift (e.g., Ivison et al. 2000b) 
and thus their bolometric luminosities often exceed 
10 12 i Q , suggesting a very high star formation rate, 
i.e., ;> 100M Q y _1 . This implies that they are mas- 
sive forming galaxies enshrouded by dusty clouds. 



(2) Lya emitters have been also found at high redshift 
z £ 2 (Hu & McMahon 1996; Cowie & Hu 1998; 
Keel et al. 1999; Steidel et al. 2000). All these are 
indeed well-defined candidates of forming galaxies 
with a high star formation rate. A remarkable class 
of such Lya emitters is called as "Lya blobs" (Stei- 
del et al. 2000) which are spatially extended ionized 
nebulae. These blobs may be interpreted as possi- 
ble evidence for the superwind activity (Taniguchi & 
Shioya 2000). 

Therefore, it seems reasonable to accept that the super- 
wind activity works at high redshift. Although it is un- 
likely that massive galaxies could be born ubiquitously in 

Hogg et al. (1998) made a comparison between various photometric redshift predictions and spectroscopic redshifts for galaxies in the 
HDF. They found that the best photometric redshift schemes predict spectroscopic redshifts with a redshift accuracy of Az < 0.1 for more than 
68% of galaxies and with Az < 0.3 for 90 - 100%. The photometric redshift for 3-379, 

%ihot — 1.72, appears significantly different from the 
spectroscopic one, z sp0 c = 2.8 (see Table 2). However, when we use a certain set of SEDs generated by using the stellar population synthesis 
code GISSEL96 (Bruzual & Chariot 1993), one obtains z phot ~ 2.2. 
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the early universe, it seems worthwhile investigating the 
gravitational instability in shocked shells driven by super- 
winds from forming galaxies located in some over-dense 
regions at high redshift; such over-dense regions beyond 
z = 3 have been found by Steidel et al. (2000) and Ivison 
et al. (2000a). 

It is also noted that we are not simply resurrecting the 
explosion scenarios proposed by Ostriker & Cowie (1981) 
and Ikeuchi (1981). What we assume here is that some 
massive galaxies were assembled from gaseous clumps at 
high redshift (i.e., the hierarchical clustering scenario) and 
then intense starbursts occurred in them, leading to the 
superwind activity in such over-dense regions. 

3.2. Model 

First we give a summary of our model. Step I: Suc- 
cessive merging of subgalactic gas clumps results in the 
formation of a galaxy with a mass of 10 11_12 M© at red- 
shift z ~ 5. Step II: Subsequently, supernova explosions 
occur inside the galaxy and then blow out as a galactic 
wind (or a superwind). This wind expands into the in- 
tergalactic space and then causes a large-scale shell with a 
radius of several hundreds of kpc. Since this radius may be 
smaller than the typical separation between galaxies, inter- 
actions of shells may also occur, resulting in the formation 
of a large-scale gaseous slab. Step III: Since the shell or 
the slab can be regarded as a gaseous sheet, filament-like 
gravitational instability is expected to occur. Step IV: 
Further gravitational instability occurs in each filament, 
leading to intense star formation along the filament. This 
is the chain galaxy phase. Step V: The filament collapses 
gravitationally into one spheroidal system like an ellipti- 
cal galaxy within one dynamical timescale of the filament 
(~ 10 8 yr). A schematic illustration of our model is shown 
in Figure 2. 

We begin our scenario with the superwind activity ex- 
pected from forming galaxies at high redshift (e.g., Ari- 
moto & Yoshii 1987). It is considered that intense star 
formation (i.e., a starburst) occurs at the epoch of galaxy 
formation in the galaxy center, producing a galactic wind 
which lasts for a characteristic time iGw(~ 0.5 Gyr for 
an elliptical with a stellar mass of 1O 11 M : Arimoto & 
Yoshii 1987). Since infalling gas is accreting onto the 
galaxy at times t > £gW) the wind interacts with this gas, 
and shocked gaseous shells form in the outer regions of the 
galaxy. If the shells are unstable gravitationally, clumps 
may be formed within them. This explosion-driven for- 
mation mechanism has been applied to the formation of 
galaxies (Ostriker & Cowie 1981; Ikeuchi 1981; Ikeuchi, 
Tomisaka, & Ostriker 1983), dwarf galaxies (Mori, Yoshii, 
& Nomoto 1999), and globular clusters (Taniguchi, Tren- 
tham, & Ikeuchi 1999). 

A strong observational constraint on such explosion- 
driven formation scenarios arises from the very small fluc- 
tuation of the cosmic microwave background (CMB) ra- 
diation observed by COBE (Fixen et al. 1996); i.e., the 
dimensionless cosmological distortion parameters are lim- 
ited to \y\ < 1.5 x 10" 5 and \n\ < 9 x 10~ 5 (95% confidence 
level) 2 . These limits lead to the fractional energy limit of 
AU/U = 4y < 6 x 10 -5 where U is the energy of the 



CMB and AJ7 is the energy converted from other forms. 
The energy of the CMB per L* galaxy is estimated to be 
- 10 63 (1 + z) ergs (Wright et al. 1994). Therefore, the 
above fractional energy limit implies that the explosion 
energy per L* galaxy must be smaller than 6 x 10 58 (1 + z) 
ergs. Since the explosion-driven formation of galaxies pro- 
posed by Ostriker & Cowie (1981) and Ikeuchi (1981) vi- 
olates this limit, they are rejected by the COBE observa- 
tions (Wright et al. 1994). 

Therefore, in our model, we take the fractional energy 
limit into account. Assuming that the superwind activity 
occurs at z = 5, we obtain the upper limit of the en- 
ergy ejected from an L* galaxy, 3.6 x 10 59 ergs. This 
energy corresponds to that released from 3.6 x 10 s su- 
pernovae because the energy released by a single super- 
nova is 10 51 ergs. Thus the total number of supernovae 
allowed per L* galaxy is Asn = 3.6 x 10 8 . Here we as- 
sume that the initial starburst occurs obeying a Salpeter 
initial mass function with the lower-mass and higher-mass 
cutoffs, Mi = 0.1 M Q and M u = 100 M , respectively. It 
is known that massive stars with a mass of to > 8M are 
progenitors of type II supernovae. In the above star for- 
mation, the number of such high mass stars per unit gas 
mass (1M©) is estimated to be ~ 0.007. Therefore, the 
mass of a galaxy which produce iVsN = 3.6 x 10 8 is esti- 
mated as M gal ~ Asn/0.007 - 5 x 10 10 M q , being roughly 
comparable to that of an L* galaxy. 

Then we consider the physical conditions of the shell 
formed by the superwind from such an L* galaxy following 
Ostriker & Cowie (1981). In the radiative cooling epoch, 
the radius of the cooled shell is 
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where E§i is the explosion energy in units of 10 61 ergs. If 
we assume that a shell is formed from this swept-out gas, 
its surface mass density is 



^cool 
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(4) 

Adopting Q = 0.01, h = 1, z = 5, and E 61 = 0.036, 
we obtain i? cool = 0.81 Mpc, M coo] = 9.1 x 10 n M Q , and 
Scooi = 2.3 x 10~ 5 g cm~ 2 for the shell formed by the 
superwind from the L* galaxy. 

Since the cooling time is estimated as ~ 2.1 Gyr, the 
shell will form at redshift z s heii ~ 1.6 given the above cos- 
mological parameters. Since the number density of galax- 
ies in the present-day universe is nc(0) ~ 0.01 Mpc -3 , 
tt-g(z) is estimated as nc(0)(l + z) 3 . Given the adopted 
cosmology, we obtain nc(l-6) ~ 0.18 Mpc -3 , giving an 



2 Recently, Battistelli, Fulcoli, & Macculi (2000) re-estimated the external calibrator (XCAL) emissivity and then obtained the following 
new upper limits; \y\ < 3.1 X 10~ 5 and //| < 4.5 X 10~ 4 (95% confidence level). These values are closer to the previous estimates by Mather 
et al. (1994). However, in our paper, we adopt the strongest constraints obtained by Fixen et al. (1996). 
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average separation, Lq = n G 1 ^ 3 (47r/3)^ 1 / 3 ~ 1.1 Mpc. 
Although this is slightly larger than i? coo i, it is suggested 
that such superwinds will overlap in some cases because 
we assume that the superwinds occur from giant galax- 
ies in over-dense regions of the universe. In either case, 
a large-scale shocked shell with a radius of ~ 1 Mpc can 
be formed in this model. The physical parameters of the 
shell are summarized in Table 3 for Q = 0.01, 0.1, and 1 
and h — 0.5, 0.75, and 1. 

Then we investigate what happens in such a shell. Since 
the radius of the shell is so large, the shell can be approx- 
imated on sufficiently small scales (i.e., ~ 10 kpc) by a 
planer, self-gravitating layer, or sheet (e.g., Elmegreen & 
Elmcgrccn 1978). It is known that gravitational insta- 
bility in an isothermal self-gravitating sheet results in the 
formation of filaments (Miyama, Narita, & Hayashi 1987a, 
1987b; Curry 2000 and references therein). According to 
Curry (2000), the critical wavelength of the filamentary 
instability is defined as 

A cr = c s [27r/(Gp shcll )] 1/2 (5) 

where c s is the sound speed, /9 s hcii is the mass density of 
a sheet and G is the gravitational constant. If we assume 
that pshdi ~ 10p cr (z) where p CI (z) = p CT (0)(l + z) 3 , we 
obtain p s heii — 3.34 x 10~ 27 g cm" 3 at z = 1.6. Since it is 
expected that the sound velocity is c s ~ 10 km s -1 (e.g., 
Ostriker & Cowie 1981), we obtain A cr ~ 54 kpc. Using 
the mass density together with the surface mass density 
of the shell, we obtain the thickness of the shell (sheet), 
AR = S coo i/p ~ 2.2 kpc. Thus an initial major-to-minor 
axial ratio is ~ 25, being consistent with the axial ratios es- 
timated for the chain galaxies (see Dalcanton & Shectman 
1996). According to Miyama et al. (1987a, 1987b), it is 
expected that this filament experiences the beads instabil- 
ity and then collapses along the major axis of the filament, 
forming an elliptical-like galaxy. Even if a cigar-like, pro- 
late galaxy could be formed from the filament, it would be 
dynamically unstable and thus evolve to a rounder ellipti- 
cal galaxy (Merritt & Hcrnquist 1991). 

The mass density of the shell leads to the dynamical 
timescale of the sheet, Td yn ~ (G/9 s hcii) ^ 2 — 2.1 x 10 9 yr. 
Since this timescale is comparable to the cooling time, it 
is considered that chain galaxies can be formed at z ~ 1.6. 
However, we note that chain galaxies could be born prefer- 
entially in higher-density regions in the shocked shell and 
thus the formation epoch may be earlier than the above 
rcdshift. In order to explain the formation of chain galaxies 
at z <~ 3, one may assume that the epoch of the superwind 
activity occurs at z > 5. 

Now let us estimate the mass of the filament consid- 
ered above. The filament could grow in mass by accu- 
mulating gas clouds within a roughly regular square re- 
gion with a length of A cr and a thickness of AR. There- 
fore, we estimate the mass of the filament, Mfy ament ~ 
A 2 r AR p ~ 3 x 10 8 M Q . Since a typical luminosity- 
to-mass ratio of O stars is ~ 10 3 , this filament becomes 
as luminous as ~ 10 11 L Q . Note that this luminosity is 
achieved if a burst of massive-star formation occurs in it 
with very high star-formation efficiency (approximately 
100%). The reason for such high efficiency is that a 
filament-like structure is inherently unstable and will dis- 
rupt on a time scale comparable to the transverse cross- 



ing time; r cross ~ 0.72 x 10 8 (d/2.4 kpc)(w/10 km s -1 ) -1 yr 
(C95). The above luminosity is roughly comparable to the 
observed luminosities of the chain galaxies although they 
show a large scatter from ~7x 10 8 L Q to — 4 x 1O 1O L 
in Ly (see Table 2); this scatter may be attributed to 
the density fluctuation in the shocked shell. After the in- 
tense star formation ceases (<~ 10 8 years after the forma- 
tion; Abraham 1997), this filament rapidly fades and then 
evolves dynamically to a dwarf elliptical-like galaxy with 
<J 1O 8 L0. This may explain why we cannot find remnants 
of the chain galaxies. 

Finally we estimate how many chain galaxies are formed 
in the shocked shell. Since the mass of the shocked shell 
is M coo i ~ 9.1 x lO n M0 and that of a chain galaxy is 
A^chain ~ 3 x 10 8 Mq, the maximum number of dwarf 
galaxies is A chain ~ 3000?y c ham,ioo where 77 chaini i o is the 
formation efficiency of chain galaxies in the shell in units 
of 100%. Since the radius of the shocked shell is i? coo i — 
0.81 Mpc, the volume number density of chain galaxies is 
estimated as n c hain ~ 1350?7 c h a in,ioo Mpc -3 . However, it is 
unlikely that galaxies were formed with such the very high 
efficiency because the chain galaxy formation could occur 
preferentially in higher-density regions in the shell. If we 
adopt r? c hain,ioo = 0.1 (i.e., 10% efficiency; see section 4.2), 
we obtain n c h a in ~ 135 Mpc -3 . 

Given ry c hain,ioo = 0.1, 300 chain galaxies could be 
formed in one shocked shell. The shocked shell radius, 
-Rcooi — 0.81 Mpc, corresponds to an angular size of 2.5 
arcmin at z ~ 1. This leads to a nominal surface num- 
ber density of chain galaxies formed in one shocked shell 
Cchain — 300/(7r2.5 2 ) ~ 15 arcmin -2 . However, since 
the lifetime of the chain galaxy phase (T c h a i n ~ 10 8 yr) 
is much shorter than the Hubble time (rnubbio ~ 10 10 
yr), the surface number density of chain galaxies made in 
one shocked shell is statistically estimated as <7 c h a i n ,stat ~ 
<7chainTchainTj^ bble ~ 0.15 arcmin -2 . As summarized in 
Table 1, the observed surface density of the chain galaxies 
is a few chain galaxies arcmin -2 . Comparing this with the 
above estimate, we suggest that ten shocked shells at dif- 
ferent redshifts (e.g., z ~ 0.5 - 3) are included in the sky 
areas of the deep survey fields. This also implies that only 
a small number of galaxies could experience the superwind 
activity leading to the formation of a shocked shell. 

Our model also predicts that several chain galaxies can 
be found at nearly the same redshift because they are as- 
sociated with one shocked shell. It is again interesting to 
note that the redshift of the Hot Dog galaxy (z = 2.803) is 
very close to the chain galaxy 3-379 in vdB96 (z = 2.775). 
Among the eight galaxies with their photometric redshifts 
(Table 2), the following two pairs have the same redshifts; 
1) Zphot = 0.92 for 2-234 and 3-621, and 2) z phot = 1.60 for 
4-241 and 3-531. Although most chain galaxies are unfor- 
tunately very faint [Tgi^AB) ;> 23], it should be possible 
to perform optical/NIR spectroscopy of the chain galaxies 
found by C95 and vdB96 in the near future. 

4. DISCUSSION 

4.1. How Can Parent Galaxies of Chains Be Seen ? 

Our model may be one of possible mechanisms to form 
chain galaxies at high redshift. Our model predict that 
there should be massive parent galaxies located within 
1 Mpc which produced the shells during a burst of star 
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formation at much higher redshift. They can be seen as 
passively-evolving galaxies at nearly the same redshift as 
that of certain chain galaxies. Using photometric redshifts 
of the HDF galaxies provided by Fernandez-Soto et al. 
(1999), we try to search for such parent galaxies. 

The most interesting chain galaxy in the vdB96 sam- 
ple is the chain galaxy 3-367 at z = 2.775 because the 
Hot Dog galaxy is also located at nearly the same red- 
shift. However, since Fernandez-Soto et al. (1999) give its 
photometric redshift z p hot = 1-72 (see Table 2), it seems 
difficult to search for its parent galaxy candidates. Then, 
as a trial, we try to search for parent galaxy candidates for 
the following two pairs of chain galaxies; 1) vdB96 2-234 
and 3-621 at z phot = 0.92, and 2 ) vdB96 4-241 and 3-531 
at Zphot = 1-60. 

Since the bin size of photometric redshift studied by 
Fernandez-Soto et al. (1999) is Sz = 0.02, we try to find 
galaxies with z p hot = 0.92 or 1.60. We have found 40 
galaxies at z = 0.92 and 35 ones at z = 1.60 in the HDF. 
However, parent galaxies should be passively-evolving ones 
and thus we reject galaxies with blue SEDs such as those 
of irregular galaxies. Then we find three possible candi- 
dates whose SEDs are relatively red at z — 0.92 and only 
one at z — 1.60. The results are summarized in Table 4. 
Among the three galaxies at z — 0.92, two galaxies (Nos. 
531 and 536) have lower luminosities than L* . Therefore, 
the remaining galaxy (No. 466) may be the most probable 
candidate although its projected separation is 0.58 Mpc 
from vdB96 3-621 and 0.73 Mpc from vdB96 2-234. As for 
the candidate galaxy at z — 1.60, its luminosity is higher 
than L* and the projected separation is 0.26 Mpc from 
vdB96 3-531 and 0.25 Mpc from vdB96 4-241. Therefore, 
this seems a promising candidate of the parent galaxy. 

Here it is again noted that the accuracy of photometric 
redshifts is not high enough to identify the parent galax- 
ies conclusively. In order to find firmer candidates, we 
will need the spectroscopic redshift for a large number of 
galaxies in the HDF. 

4.2. Connection to the Formation of Dwarf Galaxies 

As described in section 3, it is expected that the chain 
galaxies formed in our model will evolve to dwarf galaxies 
because their masses are of the order of ~ 10 8 M Q . There- 
fore, our scenario may also be related to the formation of 
dwarf galaxies. 

First, we estimate the maximum number of dwarf galax- 
ies formed in the shocked shell. Since we assume that 
dwarf galaxies are the descendents of chain galaxies, the 
maximum number is the same as that of the chain galax- 
ies; i.e., -/V(j war f = A c hain ~ 3000?7 c hain,100- 

Thus, the 

volume number density of dwarf galaxies is estimated as 

"-dwarf = "chain ~ 1350?7 c hain,100 MpC~ 3 . 

Recently, Gnedin (2000) investigated an average star 
formation law during the first 3 Gyr using the observa- 
tional data of Local Group dwarf galaxies. His analy- 
sis showed that the data are consistent with the ortho- 
dox Schmidt law with a star formation efficiency of about 
4%. Since our model is applicable to galaxies in over-dense 
regions, the formation efficiency may be higher than the 
above estimate. Therefore, it seems reasonable to assume 
that the formation of chain galaxies in the shocked shell 
is roughly 10%, as adopted in section 3. Then, the vol- 
ume number density of dwarf galaxies is down to ndwarf ~ 



135 Mpc -3 . This is significantly lower than the observed 
number density of dwarf (elliptical) galaxies in the Virgo 
cluster (i.e., an over-dense region in the local universe) is 
- 300 Mpc~ 3 (Binggeli, Sandage, & Tammann 1988). Al- 
though chain galaxies could contribute to the formation of 
dwarf galaxies, the majority of the dwarf galaxies in the 
Virgo cluster may be formed from other mechanisms; e.g., 
the gravitational collapse of less-massive gas clumps or the 
tidally-induced formation (e.g., Okazaki & Taniguchi 2000 
and references therein). 

4.3. Connection to Lyman Limit Systems and Damped 
Lyct Absorbers 

The shocked gaseous shell formed by the superwind ac- 
tivity may be important absorbing material in the inter- 
galactic space. Indeed, the surface mass density of the 
shocked shell, S coo i — 2.3 x 10~ 5 g cm~ 2 , corresponds to 
the H I column density, Ahi — 1.4 x 10 19 cm~ 2 . Since the 
temperature of IGM may be ~ 10 4 K at most, the ioniza- 
tion correction seems modest. Therefore, shocked shells 
formed by forming massive galaxies may cause some Ly- 
man limit systems (Ahi ;> 10 17 cm~ 2 : e.g., Steidel 1990) 
found in high-z quasar spectra. Furthermore, if we see 
such shells from a highly inclined angle (e.g., # v icw ;> 80°) 
or their tangential sections, the H I column density ex- 
ceeds 10 20 cm~ 2 , causing damped Lya systems (e.g., Pee- 
bles 1993; Wolfe et al. 1995). Or, more likely, the damped 
Lya systems may arise in the last stages of gravitational 
collapse of the filaments just as star formation commences. 
Small cloudlets breaking into intergalactic medium (IGM) 
may be observed as Lya forests with Ahi <i 10 15 cm~ 2 . 

The shock velocity is given by 

v cool = 200S67 5 [h 2 n (l + z) 3 ] °- U km s- 1 (6) 

(Ostriker & Cowie 1981). Adopting E 61 = 0.036, h = 1, 
tto = 0.01, and z = 1, we obtain w coo i — 130 km s _1 . 
Therefore, the absorption feature caused by one shocked 
shell may have the line width of <; 260 km s _1 . However, 
the velocity widths observed for most high-column-density 
absorption systems (determined from unsaturated low-ion 
species) are much less than 260 km s _1 . It is also noted 
that our model is applicable to the case that an absorb- 
ing galaxy is located significantly far from the quasar; i. 
e., the projected distance up to several hundreds kpc is 
acceptable. 

In our model, some giant galaxies in over-dense regions 
in the universe could experience initial starbursts and then 
show the superwind activity. Therefore, shocked shells 
may not be formed ubiquitously in high-z universe. As 
mentioned in section 3, the observed surface density of 
chain galaxies suggests that there are only several shocked 
shells at most in the redshift interval between 0.5 and 3. 
Therefore, our model predicts a few high-z Lyman limit or 
damped Lya systems per unit redshift along a ID sight- 
line." Stenglcr-Larrea et al. (1995) showed that the num- 
ber density of Lyman limit systems as a function of redshift 
is described by N(z) ~ A (l + z) 7 where A = 0.25 and 
7 = 1.5. This gives A(3) ~ 2, being consistent with our 
prediction. 

Another interesting issue related to our model is the 
metal enrichment in the IGM. One important consequence 
of our model is that such H I absorber must be polluted 
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chemically because the superwind contains a lot of heavy 
elements. It has been argued that supernova explosions 
lead to the chemical enrichment of IGM at high redshift 
(e.g., Ostriker & Gnedin 1996; Miralda-Escude & Rees 
1997; Aguirre et al. 2000). We estimate the metal enrich- 
ment due to the superwind. The mass of metal ejected 
from a star is mz — ezm* where ez is the mass fraction 
of metal with respect to the stellar mass. Thus the total 
metal enrichment due to the superwind can be estimated 
as 



AZ = A * 4 ' 2 >< ( 7 ) 

4irRt ool ARp shcn 



where we adopt to* = 10A/q as the typical mass of pro- 
genitors of type II supernovae. We then obtain AZ ~ 
4.2 x 10~ 4 if ez = 0.1. Since the metal abundance ob- 
served in Lya forests (e.g., Pettini et al. 1997; Songaila & 
Cowie 1996; Lu et al. 1998) is of the order of 0.01 Z Q where 
Zq is the solar metal abundance in mass (Z = 0.02); i.e., 
Ziqm — 2 x 10~ 4 . Therefore, the metal ejection as a result 
of the superwind from a forming massive galaxy may be 
responsible for the observed metal abundance in the IGM 
at high redshift. 

5. CONCLUDING REMARKS 

The chain galaxies found in the recent optical deep sur- 
veys have been investigated in terms of the gravitational 
instability of shocked shells formed by a superwind or by 
interaction between superwinds of forming galaxies at high 
redshift. Our new model appears consistent with the ob- 
servational properties of the chain galaxies. Since the 
chain galaxies are expected to evolve to dwarf galaxies 
with a mass of ~ 10 8 Af Q , our model explains why there is 
no local counterpart. 

We have also discussed that shocked shells may be re- 
sponsible for some Lyman limit systems and damped Lya 
systems observed in quasar spectra because their H I col- 
umn densities exceed Nm ;> 10 17 cm -2 . The metal enrich- 
ment up to Z ~ 0.01Z Q can be achieved by the pollution 
from supernova ejecta. In order to investigate this possi- 
bility further, we need very high-resolution, emission-line 
imaging of some counterparts of damped Lya systems at 
high redshift. 

The size of such shocked shells is too large (i.e., ~ 
1 Mpc) to be detected if the sky coverage of deep sur- 
veys is not so large. In order to find firm observa- 
tional evidence for such shocked shells, wide-field, deep 
surveys will be also very important in future. The 
SUPRIME-Cam (the Subaru prime focus camera; e.g., 
http: / /SubaruTclcscope.org/Introduction/instrument.htm| ) 



seems to be the most efficient tool to perform such surveys. 
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Table 1 

a summary of the surveys by c95 and vdb96 





C95 


vdB96 


Field 


SSA 13 & SSA 22 


HDF 


Area 


10.7 arcmin 2 


5.3 arcmin 2 


^Vchain(/814 < 25) a 


24 


11 


^Vchain(/814 > 25) b 


4 




0"chain(^814 < 25) c 


2.2 arcmin -2 


2.1 arcmin -2 


0"chain(-^814 > 25) d 


0.4 arcmin -2 




< a/b > c 


4.8 ±2.1 


3.1 ± 1.3 



a The observed number of chain galaxies with Jgi4 < 25. 
b The observed number of chain galaxies with igi4 > 25. 
c The surface density of chain galaxies with I S i4 < 25. 
d The surface density of chain galaxies with Jgi4 > 25. 
e The average value of apparent major-to-minor axial ratio. 



Table 2 

Observational properties of the vdB96 chain galaxies 



ID a 


ID b 


n c 


d d 


a/b e 


AB(814) f 


7 S 
^spcc 


2-phot 


M v i 


w 






(") 


(kpc) 




(mag) 






(mag) 


(£©) 


4-627 


424 


0.16 




1.12 


23.28 










3-367 


363 


0.24 




2.63 


24.43 










2-234 


1027 


0.55 


2.90 


2.63 


23.88 




0.920 


-18.5 


2.1 x 10 9 


4-375 


81 


0.37 


1.69 


3.33 


23.94 




0.640 


-17.3 


7.0 x 10 8 


4-241 


679 


0.55 


3.39 


2.63 


23.94 




1.600 


-20.5 


1.4 x 10 10 


3-243 


496 


0.45 


1.85 


2.50 


24.51 




0.520 


-16.2 


2.6 x 10 8 


3-621 


16 


0.49 


2.59 


5.26 


24.93 




0.920 


-17.4 


7.6 x 10 8 


3-379 


334 


0.34 


2.27 


2.32 


24.71 


2.775 


1.720 


-21.5 


3.5 x 10 10 


2-526 


836 


0.25 


1.63 


2.38 


25.27 




2.240 


-20.2 


1.0 x 10 10 


3-448 


276 


0.20 


1.24 


2.88 


25.35 




1.640 


-19.4 


5.1 x 10 9 


3-531 


120 


0.61 


3.75 


5.88 


24.27 




1.600 


-20.2 


1.0 x 10 10 



a ID number given in vdB96 

b ID number given in Fernandez-Soto et al. (1999) 

c Apparent size along the major axis 

d Linear size along the major axis 

c Apparent major-to-minor axial ratio 

' AB magnitude at I band 

s Spectroscopic redshift 

h Photometric redshift taken from Fernandez- Soto et al. (1999) 
'Absolute V magnitude 
3 V band luminosity 
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Table 3 
Parameters of the shell 



^0 


h 


^cool 


-f^cool 


M cool 


^cool 


AR 






(Gyr) 


(Mpc) 


(M Q ) 


(g cm" 2 ) 


(kpc) 


0.01 


1 


2.1 


0.81 


9.1 x 10 11 


2.3 x 10" 5 


2.2 


0.1 


1 


0.66 


0.32 


5.6 x 10 11 


9.0 x 10" 5 


8.7 


1.0 


1 


0.21 


0.13 


3.5 x 10 11 


3.5 x 10~ 4 


34 


0.01 


0.75 


2.8 


1.0 


1.0 x 10 12 


1.6 x 10" 5 


1.6 


0.1 


0.75 


0.89 


0.41 


6.4 x 10 11 


6.4 x 10" 5 


6.2 


1.0 


0.75 


0.28 


0.16 


3.9 x 10 11 


2.5 x 10~ 4 


24 


0.01 


0.5 


4.2 


1.4 


1.2 x 10 12 


1.0 x 10" 5 


0.99 


0.1 


0.5 


1.3 


0.56 


7.5 x 10 11 


4.0 x 10" 5 


3.8 


1.0 


0.5 


0.42 


0.22 


4.6 x 10 11 


1.5 x 10~ 4 


15 



Table 4 
Parent galaxy candidates 



ID a 


SED b 


M v c 


Lv d 




£>i e 










(mag) 


(L & ) 


(") 


(Mpc) 


(") 


(Mpc) 


z=0.92 


466 


Sbc 


-20.7 


1.67 x 10 iU 


109 


0.58 


138 


0.73 


531 


M31 


-18.8 


2.73 x 10 9 


62 


0.33 


90 


0.48 


567 


M31 


-17.8 


1.14 x 10 9 


70 


0.37 


91 


0.48 


z=1.60 


306 


Sbc 


-21.0 


2.05 x 10 iU 


43 


0.26 


41 


0.25 



a ID number given in Fernandez-Soto et al. (1999) 
b Best-fit SED type in Fernandez-Soto et al. (1999) 
c Absolute V magnitude 
d V band luminosity 

e Distance from 3-621 (z=0.92) or 3-531 (z=1.60) 
f Distance from 2-234 (z=0.92) or 4-241 (z=1.60) 
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Fig. 1. — Frequency distributions of the redshifts, sizes, and absolute magnitudes of the nine chain galaxies in vdB96. The galaxy 3-379 is 
shown by the filled square in each panel because it has the spectroscopic rcdshift. 
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Fig. 2. — A schematic illustration of the proposed formation mechanism of chain galaxies 
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